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Abstract 

We compute the 0{as) corrections to the Wilson coefficient of the kinetic operator 
in inclusive semileptonic B decays. Our analytic calculation agrees with reparame- 
terization invariance and with previous numerical results and paves the way to the 
calculation of analogous corrections to other power-suppressed coefficients. 



1 Introduction 



The precision determination of the Cabibbo-Kobayashi-Maskawa quark mixing matrix re- 
mains a central goal in the flavor physics program. A new generation of high-luminosity 
B factories is expected to start operations in a few years and in view of the improved 
experimental resolution the theoretical uncertainties should be reduced whenever possible. 
In the case of inclusive semileptonic B decays, which currently provide the most precise 
determination of \Vcb\ and \Vub\, theoretical uncertainties are already dominant, but there 
is space for improvement. 

As is well known, the theoretical foundation for our understanding of inclusive semilep- 
tonic decays B — ?■ Xciiy is an Operator Product Expansion (OPE) which ensures that 
non-perturbative effects are suppressed by at least two powers of the bottom mass m^. 
They are parameterized by a limited number of matrix elements of local operators which 
can be extracted from experimental data. The total inclusive width and the first few mo- 
ments of the kinematic distributions can be well approximated by a double series in 
and Aqco/^fe [D 12] • After extracting the most important non-perturbative parameters, 
including the heavy quark masses, from the moments, one can therefore use them in the 
OPE expression for the total semileptonic width and determine \Vcb\ from the comparison 
with the experimental rate. 

It is worth emphasizing that the information obtained from the fits to the moments of 
B — > Xc£z/ (see [21 S] for recent accounts) find other important applications. Indeed, the b 
quark mass and the OPE expectation values obtained from the moments are crucial inputs 
in the determination of \Vub\ from inclusive semileptonic charmless decays, see e.g. [5] and 
refs. therein. The heavy quark masses and the OPE parameters are also relevant for a 
precise calculation of other inclusive decays like B — )■ Xs7 [6]. 

The reliability of the inclusive method rests on our ability to control the higher order 
contributions in the OPE. If we neglect perturbative corrections, i.e. if we work at tree- 
level, we presently know the 0(l/m^) and 0(l/m^) contributions [7], while the 0(l/m^) and 
0{l/mQ) effects have been studied in [8]. Unfortunately, new non-perturbative parameters 
appear at each order in the OPE: as many as nine new expectation values appear at 
0(l/m^). As a result, only the parameters associated with the 0(l/m^'^) corrections are 
routinely fitted from experiment. In [S] the parameters associated with 0{l/mQ^) effects 
have been estimated in the ground state saturation approximation, finding a relatively 
small -1-0.4% net effect on \Vcb\. Recently, the validity of the ground state saturation has 
been investigated [9J , and it has been shown that the non-factorizable contributions can be 
in general comparable to the factorizable ones. Additional work is therefore necessary to 
assess the importance of higher order power effects. 

For what concerns the purely perturbative corrections to the free quark decay, they are 
known at O(a^) in all the relevant cases, namely for the width and the first few moments 
of the lepton energy and hadronic mass distributions. The complete 0{as) and 0{a1(3o) 
corrections have been computed some time ago (see [10] and refs. therein), while the remain- 
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ing two-loop corrections to the width and to the first few moments have been calculated in 
Refs. [m [121 [131 [S] . The theoretical uncertainty due to missing purely perturbative effects 
is now relatively small [Ti] . 

The 0(asAQ(2;£)/m^) corrections appear to be a potentially more important source of the- 
oretical uncertainty. The 0{as) corrections to the Wilson coefficient of the kinetic operator 
have been computed numerically in [15]; they can be also obtained from the parton level 
0{as) result using reparameterization invariance relations [H [161 [13 HE]- They lead to 
numerically modest 0(as/i^/m^) corrections to the width and moments, where fi^ is the 
matrix element of the kinetic operator. However, in order to have all the 0{asA.QQY)/''T^h) 
effects one should also consider the 0{as) corrections to the Wilson coefficient of the chro- 
momagnetic operator. A complete 0{asA.Q(-,Y)/''^b) calculation has been performed in the 
simpler case of inclusive radiative decays [Ldj, where the 0{as) corrections increase the 
coefficient of Hq, the matrix element of the chromomagnetic operator, by almost 20% in 
the rate. 

In this paper we present the first part of an analytic calculation of the 0{asAQQY)/ml) 
corrections. We extend the method developed in |T9] to semileptonic decays into hadronic 
final states containing a massive quark and validate it rederiving the O^asf^l/ml) correc- 
tions and reproducing the reparameterization relations. As we did in [19], we compute 
the relevant Wilson coefficients at 0{as) by Taylor expanding off-shell amputated Green 
functions around the b quark mass shell, and by matching them onto local operators in 
Heavy Quark Effective Theory (HQET). The extension to semileptonic decays implies new 
technical difficulties, because one needs to consider the mass of the final quark and the 
non-vanishing invariant mass of the lepton pair. The integrals involved are therefore less 
divergent but more complex. We identify a small number of master integrals and express 
our results in terms of the same functions appearing in the 0{as) parton calculation. In or- 
der to compute contributions to arbitrary moments, we give explicit corrections to the three 
independent structure functions, namely of the triple differential rate. The 0{as^i%/'rnl) 
calculation is under way and will be presented in a forthcoming publication. 

The outline of this paper is as follows. In section 2 we introduce our notation and review 
the known 0(l/m^) and 0{as) corrections to the triple differential rate. Section 3 is 
devoted to a description of the calculation of the 0{asf^l./ml) contributions. Our analytic 
results can be found in section 4. In section 5 we review the reparameterization invariance 
constraints and show that they are satisfied by our results. In section 6 we summarize 
and conclude. The relevant master integrals and a few technical details are given in the 
Appendix. 

2 Leading Order results 

We start recalling the tree- level results and the 0{as) corrections to the leading, free-quark 
term of the OPE. We consider the decay oi a. B meson of four- momentum pb = Mbv 
into a lepton pair with momentum q and a hadronic final state containing a charm quark 
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with momentum jj = Pb ~ The hadronic tensor W^^ which determines the hadronic 
contribution to the differential width is given by the absorptive part of a current correlator 
in the appropriate kinematic region, 

W^\Pba) = Im^ / d'xe-^''-%B\TJ^\x)jm\B), (2.1) 
uMb J 

where = c'-j^Pih is the charged weak current. The correlator is subject to an OPE in 
terms of local operators, which at the level of the differential rate takes the form of an 
expansion in inverse powers of the energy release, whose leading term corresponds to the 
decay of a free quark. 

We generally follow the notation of Ref. (TU] and express the 6-quark decay kinematics in 
terms of the dimensionless quantities 

mi mi mi 

where p = mbV is the momentum of the b quark and 

< n < u+ = (1 - f - p and < < (1 - y/pf. (2.3) 
We will also employ the energy of the hadronic system normalized to the b mass 

E = ^{l + p + u-f). (2.4) 

The case of tree-level kinematics corresponds to m = 0; we indicate the corresponding energy 
of the hadronic final state as 

Eo = lil + p-f). (2.5) 

The normalized total leptonic energy is 

% = 1 — E from which follows u = 2 {1 — Eq — %) . (2-6) 

We also introduce a threshold factor 

X = A{ql-f)=A{E'-p-u). (2.7) 

In the case of tree-level kinematics, the threshold factor becomes Aq = 4(£'q — p). It is 
convenient to introduce a short-hand notation for the square root of A: 

t = — , to = y^. (2.8) 

2E' 2Eo ^ ' 

It is customary to decompose the hadronic tensor as follows 

mb W^^Pb, q) = -Wi g^"" + W2 v^'v" + iW-^ e^'P^v^q, + W^q^'q" + {v^'q" +v''r) , (2-9) 
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where the structure functions Wi are functions of q^,qo or equivalently of q'^,u, is the 
four- velocity of the B meson, and = q'^/rrif,. 

In the hmit of massless leptons only Wi^2,3 contribute to the decay rate and one has 
dT _ Glml\V,h\^ 



dEe d(p du 



xlq^Wi 



-9{u+ — u) X 



2El - 2E,go + y 



(2.10) 



W2 + q\2E, - go) W- 



3 ( ) 



where «+ = (! — \f^Y — p represents the kinematic boundary on u, and E^ = Ec/mh is 
the normalized charged lepton energy. Thanks to the OPE, the structure functions can be 
expanded in series of and Aqcd/'"^6- There is no term linear in AQco/^b and therefore 



2ml ' 



2ml ' 



2ml ' 



2ml " 



(2.11) 



where we have neglected terms of higher order in the expansion parameters, /i^ and [j?q are 
the i?-meson matrix elements of the only gauge-invariant dimension 5 operators that can 
be formed from the h quark and gluon fields [H [2]. The leading order coefficients are given 

by 

Wf^ = ) 6iuy, t^f ) = 2Eo, = 4, = 2. 



The tree-level nonperturbative coefficients [2] read 



(2.12) 
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(7r,0) r/^N 

= wl d[u) 


(7r,0) 

wl 




Eo), w\ 


(7r,0) 
W2 


= 0, 


w\ 


(7r,0) 
W3 


= -2, 


w\ 



(1 - Eq) 
-8(1 -Eo), 

-1(1 -i^o), 



w 



w 



w- 



1 
2 

(^,2) 



3A0, 

2 \ 
3A0, 



(2.13) 



and 



iyf'°) = ^f°)5(«) + t.f^)y(^); 



(2.14) 



w 



(Gfi) 



-1(2 -5Eo) 



(G,l) 



= 0, 



w~ 



— |(-Eo + 3-Eo + l-^o) 
3-5Eo); 

-t(l + 5Eo). 



^f'^ = |(3-5Eo) 



The perturbative corrections to the free quark decay have been computed in \1Q\ and refs. 
therein. They read 



WP = wP 



Si6{u)-2{l-Eoh) 
4 



+ 



(p + u) 



+ R^ e{u) 



(2.15) 



where Si = S + Ai and 



b - Zl^o [l2fi - hfl) - I mp — hfi, 



Ai 



P J . 



A, 



4g2 



4g2 



4g2 



A, = 0, 



and the functions Ri are given in Eqs. (2.32-2.34) of Ref. [I0]|1| The integrals /i, /i,o, -^2,0, 
and /4,o are given below in Eqs. ( ]A.6|IA.7|IA.8I) and the plus distribution is defined by 



"1" 






= lim 


u 





\ne5(u) + -e(u-e)e(l-u) 
u 



(2.16) 



or equivalently by its action on a test function f{u): 

du f{u 







u 



(2.17) 



The upper limit of integration in the rhs of fl2.17p can be chosen arbitrarily, but it is 
convenient to have it larger than the physical boundary, Ref. [10] uses as upper 
limit, and the two definitions are related by the simple expression 



'1" 




"1" 


u 


+,[10] 


u 



— In u+ 5{u) — 



d{u - u+) d{l - u) 



u 



(2.18) 



3 The calculation of O^as/J^l/ml) effects 

The four diagrams in Fig. 1 are our starting point. They contribute to the weak current 
correlator of Eq. ( 12. ip and are sufficient to determine both wj^^"^ and wj^'^'^\ while additional 
diagrams with external background gluons are necessary for the determination of W^^ ' . 
The momenta of the external b quarks and W bosons are p = m^v+k and g, respectively. We 
write down the corresponding off-shell forward amplitudes, and extract the contributions 
to 14^1,2,3 by contraction with appropriate tensor projectors. We then Taylor expand around 
= 0, i.e. around the mass-shell of the h quarks, through 0{k'^). We always work in 
d = 4 — 2e dimensions and use dimensional regularization for both ultraviolet and infrared 
divergences. The result of the Taylor expansion is reduced to scalar integrals, which are in 
turn expressed in terms of 4 independent master integrals, listed in the Appendix, using 
Integration by Parts (IBP) identities [20] . 

The forward amplitudes obtained in this way correspond to quark matrix elements of 
local operators which eventually have to be evaluated in the B meson. In particular, the 

^The variables w, Ab, and r of Ref. [10^ correspond to —2Eo, A and (1 — i)/(l + 1), respectively. 
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Figure 1: One-loop diagrams contributing to the current correlator. 



fc-dependent coefficients of the master integrals should be expressed in terms of B meson 
matrix elements of local operators. To this end it is convenient to use HQET: all 0{k) and 
0{k'^) matrix elements can be expressed in terms of 

Ai = {B{v)\bSDfK\B{v)) , A2 = {B{v)%^-^G,,a>'''K\B{v)) , (3.1) 

where = d^+iggG^^T"- is the covariant derivative, G^j, the gluon field tensor, and is the 
static quark field (see |19] and refs. therein for details). While Ai and A2 are defined in the 
asymptotic HQET regime, in practical applications one deals with /x^ = — Ai + 0(l/m{,) 
and = 3A2 + 0(l/mb), defined in terms of full QCD states at finite m^. The power 
corrections to these relations are irrelevant for our calculation. In this paper we consider 
only the terms proportional to Ai. 

The next step to compute the physical structure functions Wi consists in taking the 
imaginary part of the result. There are two kinds of contributions to the imaginary part: 
the ffist comes from the imaginary part of a charm propagator, raised to power n, outside 
of the loop, via 

1 1 f— 11" 

-Imy^ Y--V^ = (3.2) 

vr y[v — qY — + tri\ nl 

the other comes from the imaginary part of a loop integral, and is related to real gluon 
emission. In the ffist case the real part of the loop integrals is multiplied by 5('u) or 
its derivatives. It follows that the real parts of the master integrals, together with their 
derivatives wrt u, are only needed at u = 0, i.e. with partonic kinematics. The derivatives 
of the master integrals wrt u can be, in turn, re-expressed in terms of the same master 
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integrals, for example using 



du 



dqodq^ d 



fixed g2 du Oq^ dq^ 



fixed 



'2 



(g2 



dqu 



(3.3) 



As a result, the contribution to the final result coming from the real parts of the loop 
integrals can be expressed in terms of a single combination of dilogarithms and a single 
logarithm, defined as I^^q and /i in Eqs. (lA.6fA.7p . The latter are the same functions that 
appear in Eq. fl2.15p . 

All the singularities are located at the threshold, u = 0. We therefore identify in the 
master integrals the terms which potentially lead to infrared divergences and employ the 
identity 



u 



-A+Be 



A-1 
p=0 



-ly 



u 



(Be)" 



ln"M 



p\ 1 + p — A + Be '■ — ' n! 

^ ^ n=0 

valid for A > 0, where the plus distributions are defined by generalizing Eq. fl2.17p . 



(3.4) 



du 



with f^\u) = '^^j^p^ ■ These terms enter the imaginary part of the master integrals and 
control the infrared divergences due to real gluon emission. 

In the calculation of 1^1,2 the problem of the ci-dimensional definition of 75 can be avoided 
by simply anticommuting it with all 7^^ matrices. In the case of the parity violating structure 
function W^, however, one needs to proceed with care and adopt a (i-dimensional definition 
for the axial- vector current. One possibility is to follow [21] and employ the replacement 



u" 







du 



m— 1 



/(«) 



p=o ^' 



(3.5) 



(3.6) 



where e^^P'^ is a strictly 4-dimensional object. This Levi-Civita tensor is multiplied by 
another antisymmetric e tensor, necessary to extract the component of the amplitude, 
and their product can be expressed as a combination of metric tensors. The latter can then 
be taken in 4 — 2e dimensions. This definition has several advantages but it violates chiral 
symmetry and its basic Ward identities. Therefore it requires the introduction of a finite 
one-loop renormalization constant Z5 



(3.7) 



of course in addition to the wave-function renormalization of the external legs. is given 
in |21] and Refs. therein for the massless case, 



a. 



vr 



(3.8) 
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and we have checked that this resuh apphes to the case of massive quarks as well. Another 
possibility, which leads to the same result without an extra finite renormalization, is to 
anticommute 75 to the extreme left of the Dirac string in all diagrams and then to replace 
it by its four-dimensional definition. 

After combining all the diagrams, the infrared divergences are cancelled and the ultraviolet 
divergences are removed by the b quark wave function renormalization and the charm mass 
renormalization. The charm mass renormalization also removes all the 6"'{u) terms. No 
renormalization is necessary in the effective theory p9] . 



4 The 0{asi^l/ml) results 



We now report our results for the 0{as) corrections to the Wilson coefficients of the kinetic 
operator. The most singular part of the WiS has a universal structure 



R - 



S, + 3(1 - Eo/1,0) S"{u) - 4 (1 - Eo/i , 



(4.1) 



and the complete results are 
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5{u), (4.2) 



W. 



45(2,.)- 16(1 -Eo/1,0) (1-^0 



e{u) 



(4.3) 
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2B. 



i3,n) 
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"3 
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(4.4) 



Ao 



+ ^ + lj(l-£;o/i,o) + (i-£;o)/i,o 



^o(^o-p) 



5{u). 



In the above expressions the coefficients of the derivatives of S (u) have been reduced using 
integration by parts: for instance 



f{u) 5"{u) = m - 2/'(0) 5\u) + r (0) 5{u). 



(4.5) 



The coefficients of the plus distributions can be similarly reduced by Taylor expanding them 
around it = 0, for example: 



fin) 



/(O) 



+ /'(o) 
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/(«)-/(0)-«/'(0) 
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where the last term is regular for m — )■ 0. 
The functions B!f ^ are related to real gluon emission and are given by 
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4:{l-E){8iE-E^) -u{l-9E + 2u)) / Ix AE^u 2 {AE"^ + 3Eu + u) 

3A^^ \^~e)^ ~ 3? 

- 20^2 + 8z(z + 1) 16E(E - z) 7 - 3E - 2u 7E-2u-l 

H ■ — ■ — - H H \ , (4.9 

3uz'^ Su'^z Sz"^ 3Ez ' 

with z = u + p. Of course, the functions Rl are regular in the hmit m — )■ 0. 

The above resuhs can be used in Eq. (12. lip to compute the 0(as/i^/m^) corrections to 
the total rate and to the moments of various differential distributions. The phase space 
integration is rather delicate because of strong cancelations between different singular terms, 
especially in the presence of a cut on the lepton energy. We have compared the numerical 
results with those in Tables 1 and 2 of Ref. [15] and found agreement in all cases. In 
principle, it is also possible to take the limit p — )■ and obtain analytic results for the 
B — !■ Xuity decay. 



5 Reparameterization Invariance relations 



Reparameterization Invariance (RI) [Ml [IT] connects different orders in the heavy quark 
expansion. In particular, as we have mentioned in the Introduction, RI links the coefficient 
of the kinetic operator to the coefficient of the leading, dimension 3 operator. In the total 
rate this corresponds to a rescaling factor 1 — yU^/2m^ on the leading power result, which is 
nothing but the relativistic dilation factor of the lifetime of a moving quark and applies at 
any order in perturbation theory. The relations for differential distributions and moments 
tend to be more intricate, see [15], especially in the presence of experimental cuts. 

Recently Manohar has derived elegant RI relations [18] that apply directly at the level 
of the structure functions Wi. They are also valid to all orders in perturbation theory and 
give the a" coefficient of p^ in terms of the leading a" coefficient and its derivatives: 
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To verify these relations from Eq. (12.151) we need the first two derivatives of the plus 
distribution of Eq. (I2.17p . They can be re-expressed in terms of the higher order plus 
distributions introduced in Eq. (13. 5 p and of delta functions: 
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where we have neglected terms that do not contribute upon integration in the physical 
range fl2.3p . The coefficients wj^'^'^^ obtained from Eq. fl2.15p using the RI relations agree 
with the results given in the previous Section. Using Eqs. (15. ip one can also verify the 
relations between moments with and without cuts given in jTSj . 



6 Summary 

We have presented an analytic calculation of the 0{as) corrections to the Wilson coefficient 
of the kinetic operator in inclusive B — )■ Xcii' semileptonic decays, following and extending 
the method developed in Ref. ^9]. We have conffimed the numerical results presented in 
[T5] and reproduced the RI relations given by Manohar [18]. We have provided several 
details of the technique we have used; in particular, the Appendix contains all the master 
integrals. 

The calculation represents the ffist part of a complete study of the perturbative corrections 
to the coefficients of the power suppressed dimension 5 and 6 operators, and has offered us 
the opportunity to perform various checks. Our technique is currently employed to compute 
the Oi^asii'c/ml) corrections, as done already in the case of — )■ in |T9]. 



Acknowledgements 

We thank Lorenzo Magnea for useful discussions. PG is grateful to Thomas Manuel and his 
group at the University of Siegen for the warm hospitality in June 2012 and to the CERN 
Theory Unit for hospitality during the final stage of this work. Work partly supported by 
MIUR under contract 2010YJ2NYW_006 and by the DFG through the SFB/TR 9 grant 
"Computational Particle Physics" . 



Appendix 

In this Appendix we list the master integrals relevant to our calculation. Whenever ap- 
propriate, they are expanded in e up to the order which is necessary in our calculation. 
As explained above, we need both the imaginary and real parts of the master integrals. 
However, since their real parts always appear multiplied by 5{u) or its derivatives, the real 
parts and their derivatives wrt u are only necessary at m = 0, i.e. with partonic kinematics. 
We introduce 
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Let us also employ z = {p — qY /ml = it + p. The massive tadpole integral is always real 
and is given by 



^{0,1,0) - ^6 —2 



- + l + e(l + -)+0(e 



(A.2) 



with /(o,o,i) given by the same expression with f-)- rric. We have two distinct two-point 
functions. The first one develops an imaginary part for {p — qf > m^: 
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ml 
z + p 
2z 



2 \ e /•! 



dx{l — x) 



1 X — ir] 

P 



(A.3) 



^.-1^1-2. 



- + 2 - — -1 



2e 



1 u\nz + 2p\np 

e z + p 



+ 



0(6)}. 



The value and derivatives wrt u of /(i,o,i) at it = can be readily obtained from the above 
integral representation: 



dl. 



(1,0,1) 



du '"=0 



^1 p- 

mi 



+ 2 + 0(e) 



du^ l"=o 

C?^-^(1,0,1) I _ 



-3 



mi / 



m 



+ 1 + 0(e) 



i + i + 0(e) 



The second 2-point function is always real in the kinematic domain we are interested in. 
It can be directly expanded in e: 



/(o,i,i) = e'"T(e) f 4) [^P + (1 - - ^f)] 



(A.4) 



p^ 



mt 



- + 2 + ^^^1 In p + In + O e 

e 2g^ 1 — h 



The only three-point function is /(i,i,i); we reduce it in the following way 



Ai,i,i) 



— e 



7b e 



mt \ m 



/ dxdyx '[xx- {l-y){u + iv)] (A.5) 
V"^6/ Jo 



I f p 



2 \ e /.I 



^6 \^6/ Jo X 



M h - In ■ 



eV2e 2 {1 - 
-^ + lnx-^lnx + 0(e) 



+ 0(6) 
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where x = + 1/(1 — — z) + q^y"^, X = p + 1/(1 — (f — p) + q^V^ ■ 

We need the real and imaginary parts of /(i,i,i), as well as its derivatives wrt m at -u = 0. 
The latter can be computed from the Feynman parameter integral or can be expressed in 
terms of the master integrals at m = by applying the derivative wrt u at the integrand 
level in (lA.ip . see Eq. (13. 3p . We obtain, up to 0(e) terms, 



M=o 2 mi 



-/i,o - I 



4,0 



du 



^(1,1,1) 



1 



n=o Aom^ 
2 



(p - ^o)(^ - 2) + p(l - Eo)(^ - /4,o) + ^0 Inp 



u=0 



6a — /3/p^ /a 



4e 



An 



/3 



4p2 



7 1 , l-g^(l + 2p/i,o) 

— a 

2 2 2p 

where /3 = -8E^{Eo + p) + 4p(5Eo - p), a = 2(1 - Eq), and 



i-i 



1 1 
^1,0 = / dy—— 

lo x{y) 



In 



^^X(l/) ' 

X{y) to Eq 
+ lni±^lnEo(l-to)^(l+to)UlnEo(l+to)ln j~^°i|~^°| 



Li2 (ai)-Li2 (02) 



and 



ai 



2to/(l + to 



l-Eo(l-to) 
For the imaginary part we also need the integral 



02 



2to Eq 



l-Eo(l-to)' 



and its first derivatives at m = 0: 



'^^^ \n[x{y)/{l-yn 
x{y) 

= h,o ~ 2/2,0; 

M=0 



dl^ 

du 



- 2(1 - Eq) (p-^lnp 

— ^ l-'4,0 ~ ^-'1,0 ~ ^-'2,0j + ^ 7 ) 

n=0 Ao Aq P 



(A.6) 



(A.7) 



dv? 
where 



(3g2 + 2/;2t2) (j^^o _ 3/^_„ _ 2/2,0) + q\hfl - Eo/p) (5 - Qap" 



u=0 



(lnp-1) 



-^2,0 



^ ^^ ln(l - y) _ Li2(l - ^0(1 + to)) - Li2(l - ^0(1 - to)) 
X(?/) 



(Ai 
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